Engineering atomic quantum reservoirs for photons 
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We present protocols for creating entangled states of two modes of the electromagnetic field, 
by using a beam of atoms crossing microwave resonators. The atoms are driven by a transverse, 
classical field and pump correlated photons into (i) two modes of a cavity and (ii) the modes of 
two distant cavities. The protocols are based on a stochastic dynamics, characterized by random 
arrival times of the atoms and by random interaction times between atoms and cavity modes. The 
resulting effective model yields a master equation, whose steady state is an entangled state of the 
cavity modes. In this respect, the atoms act like a quantum reservoir, pulling the cavity modes into 
an entangled, Einstein-Podolski-Rosen (EPR) state, whose degree of entanglement is controlled by 
the intensity and the frequency of the transverse field. This scheme is robust against stochastic 
fluctuations in the atomic beam, and it does not require atomic detection nor velocity selection. 



I. INTRODUCTION 

The quest for quantum control of mesoscopic systems 
imposes the development of novel strategies, which go 
beyond the implementation of fully coherent Hamilto- 
nian dynamics. The main objective is to combat the 
effects of noise, which give rise to dissipation and de- 
coherence. Protocols based on quantum error correc- 
tion m [5], quantum feedback [21 H], dynamical decou- 
pling 121 E] , and decoherence-free subspaces [THS] aim at 
minimizing the detrimental effects of coupling to an ex- 
ternal environment. In this context, quantum reservoir 
engineering makes use of noise statistics as a resource for 
implementing robust quantum dynamics. The basic idea 
is to implement a stochastic dynamics whose stationary 
state is a nonclassical state. This is achieved by manip- 
ulating the coupling to the reservoir, whose properties 
are known only through the statistical averages [lOj. A 
prominent example of quantum reservoir engineering is 
laser cooling, which achieves low temperatures of single 
atoms or ions by tailoring the scattering cross section, 
such that in average the scattered photon carries away 
mechanical energy of the atomic center of mass [11] [12] . 
More recently, this concept has been generalized in var- 
ious directions, for the generation of nontrivial many- 
body states and nonequilibrium quantum phases |13j . 
and for the implementation of dissipation-driven quan- 
tum computation [14] . 

The idea of using an atomic beam as a reservoir for 
a field in a cavity actually goes back to quantum laser 
theory, see for instance Ref. [TS] , where a thermal atomic 
beam crossing the laser cavity acts as a thermal bath for 
the laser field. Actual implementations of this idea have 
been made in experiments on cavity quantum electrody- 
namics, where a beam of atoms in the lower state of a 
two-level space is used to generate a vacuum state of the 
field in a microwave cavity, by absorbing thermal photons 
initially present in the cavity [rBHTO] . 



Recently, we proposed a method for preparing quan- 
tum states of the electromagnetic field based on quantum 
reservoir engineering [20| . This method is implemented 
in a typical setup of microwave Cavity Quantum Electro- 
dynamics (CQED) as in Fig. [ija), where the resonator is 
pumped by a beam of atoms with random arrival times, 
and needs neither atomic detection, nor detailed control 
of the sequence of atoms. We showed that, by suit- 
ably preparing the initial state of the incoming atoms, 
two-mode squeezing, i.e., Einstein-Podolski-Rosen corre- 
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FIG. 1: Setup of the system for creation of EPR state of 
a) two modes of the same cavity and b) two spatially sep- 
arated cavities. In both cases atoms from a beam are first 
prepared in a coherent superposition of two Rydberg states 
\g) and |e) by a combination of laser and microwave fields. 
The atoms have random arrival times, and a low pumping 
rate warrants that at most one atom is inside the resonator 
at a time |22l I23| . While in the cavities, the dipole transi- 
tion \g) — >■ |e) is saturated by a transverse microwave field, 
thereby pumping on resonance two nondegenerate modes of 
either one or two cavites, which are led asymptotically to a 
two-mode squeezed state. 
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lations [5T] , is created between the cavity modes at steady 
state. In this paper we provide a detailed discussion of 
the proposal and analyze its robustness. We show, more- 
over, that further nonclassical states of the electromag- 
netic field can be realized by tuning different parameters 
of the setup in Fig. [ija), and which are stationary states 
of the interaction with the atomic beam. 

We also propose a scheme for entangling two distant 
resonators in a setup like the one sketched in Fig. [ijb). 
Our procedure extends to the preparation of nonlocal 
EPR states previous ideas regarding the production of 
macroscopically separated fields [24H27] . 

This article is organized as follows. In Sec. |TT] we re- 
view the basic concepts at the basis of quantum reservoir 
engineering in microwave cavity QED. In Sec. |III| we dis- 
cuss the specific scheme for creating EPR states of the 
modes of a resonator. In Sec. IIVI a method for entan- 
gling the modes of two distant cavities is proposed. The 
conclusions are presented in Sec. |V| 



II. AN ATOMIC RESERVOIR FOR 
MICROWAVE PHOTONS 

A typical setup of microwave cavity QED is sketched in 
Fig. [ija). Preparation and monitoring of the cavity field 
is achieved by interaction with atomic beams, whose in- 
ternal state is prepared in a circular Rydberg state with 
large principal number, and possessing a dipolar tran- 
sition which couples resonantly with the cavity mode. 
While the internal state of the atoms and their velocity 
can be prepared with high precision, the atomic arrival 
time is a stochastic variable which is known only prob- 
abilistically, according to a Poissonian distribution. The 
interaction of the cavity field with individual atoms is 
hence warranted by setting very low rates, such that the 
probability that two atoms are found inside the resonator 
is negligible f22l [23] . The interaction between atoms and 
cavity mode is coherent to a large extent: The atomic 
circular Rydberg states are typically stable over the in- 
teraction time, and high-finesse resonators in microwave 
cavity QED can reach very long photon-storage times, 
which in state-of-the-art experiments can reach the or- 
der of fractions of seconds [231 [2H] • These properties al- 
lowed for the realization of milestone experiments. Some 
paradigmatic ones are the preparation and measurement 
of nonclassical states of the microwave field [2H [53] , the 
experimental characterization of loss of coherence of the 
quantum field |29j and of the transition from quantum 
to classical dynamics |30| . and the quantum nondcnioli- 
tion measurement of the number of photons of the cavity 

field nil m]. 

Most of these formidable results were obtained imple- 
menting Hamiltonian dynamics, where the interaction 
between atoms and cavity photons is essentially disper- 
sive. On the contrary, the realization of nonclassical 
states of the cavity field in the dissipative regime, where 
atoms and photons exchange energy, is based on prop- 



erly tailored dissipative dynamics, where the atoms act 
as reservoir of the photonic field. As mentioned before, a 
simple example is the preparation of the cavity mode in 
the vacuum state [TBHTO] . Creation of other Fock states 
of the cavity field has been achieved in milestones ex- 
periments made in the strong-coupling regime, when the 
resonator field saturates the atomic transition [351 [33] . 
In this case, by accurately selecting the atomic velocity 
so that the interaction time of each atom is fixed, the 
nonlinear dynamics of atom and cavity field possesses 
several fixed points, so-called trapping states [31], which 
approach Fock states of the cavity field in the limit of 
negligible dissipation. 

Several theoretical proposals have been discussed in 
the literature, which provide schemes for the preparation 
of an arbitrary single-mode quantum state of the electro- 
magnetic field in a resonator, involving resonant inter- 
action with a well-controlled sequence of atoms, without 
the need of atomic detection [35H37] . This latter require- 
ment is indeed important, since present experiments lack 
high-efficiency detectors. 

In this section we briefly review the basic properties of 
atom-photon interactions in microwave CQED from the 
point of view of quantum reservoir engineering, hence set- 
ting the ground for the proposals for establishing EPR- 
correlations presented in Sec. |III| and Sec. |IV| of this ar- 
ticle. For a comprehensive review the reader is referred 
to Refs. [TS1[11[13[3H] 



A. Jaynes-Cummings Hamiltonian 

The stability of the atomic states and of the resonator 
mode during the interaction time justifies the use of 
Hamiltonian dynamics for the interaction between a sin- 
gle atom and the cavity mode. The Hamiltonian govern- 
ing the dynamics of a single dipole and the cavity mode 
is well described by the Jaynes-Cummings Hamiltonian, 
H — Hq + Hjc- Here Hq = huoO'^d' + hjJcO^a gives the 
unperturbed evolution of a dipole at frequency and 
a harmonic oscillator, the cavity mode, at frequency Wc, 
while the interaction between the dipole and the electric 
field of the cavity reads 

-ffjc = %a^a- + H.c, (1) 

with g the coupling strength. Here, a, a) are the an- 
nihilation, creation operator of a cavity photon, while 
(7 = |5)(e|, ct'^ = \e){g\ describe the dipole lowering and 
rising operators, with \g) and |e) the ground and excited 
states, respectively, of the dipolar transition. 

For wq = Wc, then [Hq,Hjq\ = and the evolution 
operator can be written as 

ij{t) ^ exp (^-iHt/hj ^ exp (^-iHot/hj exp (^-iHjct/hj 

(2) 
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whereby 



/%) = cos(0Vata)|g) + 



. sm 



«|e>(3) 




|e) = cos(0v (20^)16) + ia 



and (j) = .9^ is the Rabi angle. These equations show 
exphcitly the periodic exchange of energy between dipole 
and field when they are resonantly coupled. 

In the rest of this section we will assume that cavity 
mode and atomic dipole are resonant, iOc = wq. 



B. Interaction with a beam of atoms 

The interaction of the cavity mode with an atomic 
beam gives rise to a stochastic evolution, which is mainly 
due to the statistical knowledge of the number of atoms 
which have interacted with the cavity mode. It is then 
appropriate to use a density-matrix formalism for the 
cavity-mode state. We denote by p the cavity-mode den- 
sity matrix. The time evolution of the cavity field is 
characterized by: (i) the typical time scale which sepa- 
rates the arrival of two atoms, which is here given by the 
arrival rates r and r^, for the atoms prepared in state \g) 
and |e), respectively, (ii) the interaction time r between 
atom and resonator, which is determined by the atomic 
velocity and which follows a normal distribution V{t) 
(which we do not specify for the moment). Typically, 
rr, TeT <C 1 and one can study the field evolution on a 
coarse-grained time scale Ai, such that At ^ f, with 
f the average interaction time. For rAi,reAt ^ 1 such 
that we can assume that there is at most one atom inside 
of the cavity, we can determine the density operator at 
the time t + At given /5(t), according to the formula 

p{t + At) = p{t){l - rAt - rUAt) 

+rAt / d<j>p{(l)) 
Jq 



cos(0\/ a^a)p{t) cos(0\/ a^a) 

+71 cos{(j)\/ aat )p{t) cos{(j)\/ aat ) 
. sinf^V dd^ ) sin((/i\/ dd^ ) 

— T^^^W — 7^^' 

vaaT Vaa' 



-Ud 



(5) 



where we used — TZr and we wrote the distribution 
P(r) in terms of the distribution p{(j>) of the Rabi angle 
(j) = gr. Operator Ap in Eq. ^ is the differential change 
of the field state. The master equation is given by the 
differential equation dpt/dt which is found from equa- 
tion Ap/At. This equation has the Lindblad form [55] . 



as one can verify, but it has no trivial solution. Since the 
atoms are initially uncorrelated with the cavity mode, 
the inhomogeneous term of the Zwanzig-Nakajima mas- 
ter equation disappears [1^. Moreover, the Markov ap- 
proximation is valid in presence of a single cavity: the 
atoms exiting the resonators do not interact with it any 
longer, leading to no memory effects. Indeed, in Eq. ([5]) 
the density matrix at time t + At depends only on p(T). 
We remark that, if the atoms exiting one resonator then 
interact with a second physical system, then correlations 
mediated by the atoms must be taken into account. 

Below we discuss the master equation for the cavity 
mode in two specific limits: the weak-coupling regime, 
when the coupling of a single atom with the cavity mode 
is a small perturbation of the cavity state, i.e., 0^1, and 
the strong-coupling regime, when a single atom perturbs 
significantly the cavity state, and > 1. From now on 
we denote by pt the density matrix of the cavity field at 
time t after the coarse-grained time averaging. 



C. The weak-coupHng Hmit 

The weak-coupling limit corresponds here to the 
regime in which the mean interaction time f fulfills the re- 
lation gfy/n <C 1, where n is any relevant photon number 
(such that the corresponding populations and coherences 
are not negligibly small) and the width of the distribu- 
tion for T is assumed to be small compared to the average 
value f. In this limit the field operators in Eq. ([s]) can be 
expanded in powers of 0, and the dynamics of the cavity 
density matrix pt is governed by the master equation |17j 

^ = -iujc[d'^ d, p] - ^ {d^dpt - 2dptd^ + ptd^ d) 



7e 



[dd^ Pt ~ 2d'^ptd + ptdd^) , 



(6) 



where 7 = rg^f and 7e = TZ^f [H]. This equation de- 
scribes the incoherent energy exchange between the cav- 
ity field and an external reservoir, with loss and pump 
rates 7 and 7e, respectively. 11 TZ < 1, the resonator 
thermalizes with an effective reservoir at temperature 



T 



1 hujc 

Kb I In 72.1 



The steady state of the resonator is hence a thermal state, 
whose temperature can be controlled by adjusting the 
parameter TZ, giving the average rate of atoms prepared 
in the excited states over the ones prepared in the ground 
state. 

The steady state can be a pure state, the vacuum state, 
by preparing the atoms exclusively in the ground state \g) 
(which corresponds to setting 7?. = in Eq. Q): In this 
case the atoms absorb in average energy from the cavity 
mode until it reaches the vacuum state |0). This identifies 
a simple strategy which allows one to prepare the cavity 
in the vacuum state as a result of the interaction with a 



4 



beam of atoms, of which one controls only the initial in- 
ternal state and the arrival rate. This procedure does not 
require atomic detection nor control of the atomic veloc- 
ity, but only the mean value and the variance of the veloc- 
ity distribution, so to warrant the weak-coupling regime. 
Following the line of reasoning presented in Ref . [13] , this 
strategy could have also been identified on the basis of 
the observation that the state jg, 0) is an eigenstate of 
Hamiltonian H such that ifjc|(7,0) — 0, and it is the 
unique dark state of this dynamics. On this basis, one 
can construct a master equation which has as steady state 
\g, 0), and which has the form given in Eq. ([6]) for 7?. = 0. 

A useful benchmark for the quantum state preparation 
is given by the fidelity J^{t) for preparing the system in 
the vacuum state at time t since the beginning of the 
experiment. The fidelity T{t) = (0|p(t)|0) corresponds 
to the population of the vacuum state at time t. The 
solution p{t) can be exactly evaluated using the damping 
basis 021, and reads 



m = 



^e--*p(°)«(°) 

n=0 



(7) 



where pn^ are the right eigenvector of the Liouvillean de- 
fined in Eq. ^ corresponding to the eigenvalue 77.7, and 

an'' is the n-th moment of the expansion in the number 
operator, taken over the initial state of the cavity field. 
Using the explicit form one finds [T7] 

^W = l-(n)oe-* + Ml^e-^* + ... (8) 

where {f{a^a))o = Tv{f{a^ a)p{Q)}. When the initial 
state is thermal, with average photon number n = 
(a^ap(O)), then expression (l8| takes the compact form 



-Fwc(i) = Y 



1 



1 



-It 



1 



ne 



-<t>^rt 



(9) 



and in this case one sees that the time scale for reaching 
the ground state is determined by the damping rate 7 
and by the initial occupation number n. In other words, 
the achieved fidelity is given by initial occupation number 
n, the average Rabi angle 0, and the average number of 
atoms N = rt which have crossed the cavity. 



D. The strong-coupling limit 

The strong-coupling limit is characterized by > 1. 
In this regime, one can see from Eqs. Q and Q that, 
for a given interaction time tq = Tr/g^mo (and any in- 
teger multiple of tq), with mg > integer number, there 
is no net exchange of energy between cavity field and 
dipole when the initial state of the system is a coherent 
(or incoherent) superposition of the states \g, mo) and 
|e, mo — 1). In other words, these states are fixed points 
of the resonator dynamics, in the absence of dissipation. 



They are however not unique: Indeed, the subspace of 
the fixed points of the evolution operator for a fixed in- 
teraction time To has infinite dimension, being at least 
composed by all states \g, m), |e, m — 1), with m = ^^mo 
and £ = 1,2,.... These states have been denoted in the 
literature as trapping states [32l [34] . 

The theory of trapping states has been reported in |34j . 
These properties have been used in milestone experi- 
ments |32j in order to generate photon number states 
of the cavity field. The effect of noise on this dynamics 
have been theoretically analyzed in [?31|33]. Using simple 
considerations we now discuss how the vacuum state of 
the electromagnetic field can be the unique asymptotic 
state of the dynamics, and determine the corresponding 
preparation fidelity. For this purpose, we assume that all 
atoms are initially prepared in their ground state {TZ — 0) 
and that their velocity (interaction time) is distributed 
according to a function p{4') with finite width. The time- 
evolution of the diagonal elements of the field density ma- 
trix, {n\p{t)\n) = Cn{t) (with J2n Cn{t) = 1 as Tr{p} = 1) 
is found from Eq. ([5]), and is given by the set of coupled 
differential equations 



Cn{t) = 

with i?o = and 



Bn — 



Their solution reads 



i(i) -t- rBn+lCn+l{t) 



co{t) = co(0)+rBi / dTCi(r) 



(10) 



(11) 



(12) 



c„(0 = e-^"'-*c„(0)+ri?„+i / dTe-^"'-(*-^)c„+i(r), 

Jo 

showing that the vacuum state is always a trivial station- 
ary solution. Other stationary solutions can be found if 
there exist n > for which the coefficients Bn = 0. 

In order to study the behaviour of the coefficients i?„, 
let us assume that the velocities follow a Gaussian distri- 
bution, such that p{4>) = exp {—{(f> — 0o)^/2cr^) /V^tw^, 
where 4>o is the center and a the width. The coefficients 
B„ then read 



B„ 



1 



e-^'"" cos(20oA 



(13) 



such that for cr > only Bq vanishes. In this limit, hence, 
the only stationary state is the ground state. Neverthe- 
less, for sufficiently small values of the width a there 
exist coefficients Bn, with n > 0, whose value can be 
very small, so that the occupation of the oscillator ground 
state may converge very slowly toward unity as a function 
of time. A limiting case is found when ct — > 0, such that 
p{(l>) — S{4> — (j)o). This is the situation in which trapping 
states may exist. For instance, if we choose (f>o = 7r/2, 
then the coefficients Bq and -B4„2 (for n = 1, 2, . . .) van- 
ish, indicating that the Fock states |0) and the states of 
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the subset {|4n^)} are fixed points of the dynamics, i.e., 
trapping states. The initial state of the cavity determines 
in which of these trapping states the cavity will be found. 
In general the final state is a statistical mixture of these 
states. 

The other limiting case is found for large cr, i.e. a broad 
distribution for the atom's velocity. For tr > tt, each of 
the coefficients samples approximately equally over the 
period of the sinusoidal function in Eq. (Ill, avoiding 
trapping states, and we can therefore set Bn — \. For 
the fidelity in this approximation we find 



J-.e(t) = l-e-'--'/2^c„(0)^ 



{rt/2Y 



(14) 



Tl = 



where we used the relation T"e ^'^At 



1 



Dt X^n (Dt)" 



En 
m=0 



For an initial ther- 



mal state with c„(0) — (1 — , where /i is given 

by the initial average number of thermal photons via 
{a^ a}) = Y^-T, this simplifies to 



2^-rt(l-A'')/2 



(15) 



and one sees that in this limit the fidelity of the ground 
state is determined only by the initial occupation number 
n (parametrized by ^) and by the number of atoms rt 
which have crossed the cavity. 



III. GENERATING EPR ENTANGLED 
RADIATION WITH AN ATOMIC RESERVOIR 

In the following we discuss in detail and extend a pro- 
posal for generating EPR-entangled states of two cavity 
modes by resonant interaction with a beam of atoms, 
which was first presented in Ref. [20 . The setup is 
sketched in Fig. [T](a), and is composed by a resonator, 
with two modes in which a single atom can emit photons. 
The atoms are part of a beam crossing the resonator, and 
are driven by a classical microwave field while they inter- 
act with the cavity modes. The interaction is shaped in 
such a way that the atoms emit correlated photons into 
the cavity modes. The underlying mechanism is four- 
wave mixing, where emission into the cavity modes is 
enhanced by resonant coupling with the Rabi sidebands 
induced by a classical field that saturates the atomic 
transition [15] , with the creation of EPR-correlations be- 
ing enforced by the initial quantum state of the injected 
atoms. Contrary to typical setups based on optical para- 
metric amplifiers, the atoms pump the resonator through 
resonant single-photon processes. 

Before discussing details of the mechanism, let us make 
some general remarks on the main ideas regarding the 
generation of EPR states through reservoir engineering. 
We first label by 1 and 2 the two cavity modes, such 
that ojj is the frequency and a^, a| are the annihilation 
and creation operators of an energy quantum hujj of the 



corresponding mode (j = 1,2). Our goal is to generate a 
dynamics described by a master equation whose steady 
state is a two-mode squeezed state, pst — '^ith 



|^)=5t(O|0,0). 



(16) 



Here, |0, 0) is the vacuum state of both cavity modes and 



S'(^) = exp (Caia2 - (,a\a\ 



(17) 



is the two-mode squeezing operator, with ^ a complex pa- 
rameter. It can be verified that state is the vacuum 
state of harmonic oscillators, whose annihilation opera- 
tors hi and 62 are related to the cavity modes operators 
by the relation 



6i = 5t(e)ai5(0, 
h2 = S\i)a2S{0, 



(18a) 
(18b) 



and analogously for the creation operators. We refer to 
modes 61 and 62 as Bogoliubov modes, in analogy with 
the Bogoliubov transformation used in solid-state physics 
j46j. State p^t can be hence the dark state of an inter- 
action Hamiltonian Hi, which has, say, the form of the 
Jaynes-Cummings Hamiltonian -ffjc in Eq. ([l]) but with 
the operator a (a^) replaced by the operator bi {b\). We 
note that it is also the dark state of the interaction Hamil- 
tonian H2, which has the same form as Hi but now with 
&2 and b\ in place of bi and b\. In particular, pst is si- 
multaneously dark state of both interactions, and it is 
unique. Hence, a dynamics can be constructed, which 
has as unique steady state pst, by implementing sequen- 
tially two interactions which effectively damp oscillators 
hi and 62, respectively. In the following we will show how 
to engineer such dynamics. 



A. Engineering the coupling to the reservoir 

The Hamiltonian of driving field, atom, and cavity 
modes, in the reference frame rotating at the frequency 
ujl of the classical field, has the form 



= 1,2 j = l,2 



cra\] (19) 



where gj are the coupling constants between the two-level 
atom and each cavity mode, detunings 5j 



A = — ujQ, and 



(20) 



describes the coupling between dipole and classical field, 
with strength il, see Fig. [2] 

Let the coupling to the classical field be much stronger 
than the couphng to the cavity modes, ^ \g\\, it is 
then convenient to express Hamiltonian ( 19 ) in the basis 
of eigenstates |±) of Hq, with 



?io |±) = -fi(A T rf)/2 |±) , 



(21) 
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and 



d = v/A2 + 4rj2 



(22) 



The states |±) are the semiclassical dressed states, and 
read 



1+) = smd\g) + cos6' |e) 
|— ) = cos9\g) — sm9\e) 



with 



tan 9 



2\n\ 

d-A 



(23a) 
(23b) 



(24) 



The corresponding energy levels are shown in Fig. [2] 
We introduce the raising and lowering operators for the 
dressed states basis, 

^+ = l+>(-|, ^- = |-)(+|, (25) 

with TTz — 1+) (+1 — |— ) (— I- Using the semiclassi- 
cal dressed-state basis and operator notation we rewrite 
Eq. Iigl as H^^ = Ho + ffint, with 



Ho hdTr^/2 - ^ V Sxd\dx 



A 



and 



Hint ~ 



X 



T^z [ax + a^^ 



cos f sm ( 



(26) 



TT d\ + TT^at I sin 



If 1.9a I ^ we can choose which processes are resonant, 
and thus relevant for the dynamics, by changing the 
values of 5i,(52, and d. One is thus able to generate 
a diversity of dynamical processes, some of which are 
discussed in Refs. [201 ITf] . 

Let us set 5i — d and 82 = —d, as shown in Fig. [2] If 



\g\\ <^ d we obtain from Eq. (26) the effective Hamilto- 
nian H^^ w Hctf = Hq + Hi^t , with 



Hint — f^g (0,2 '^'^^^ ' 



ai sm 



H.c. 



(27) 



and Hq — hd{nz/2 — a\di + a2'^2)- We have assumed 
g gi = g2- The processes described by Eq. (27) are 
indicated by the arrows in Fig. [2] In the basis of the b 

operators, we have Hq — hd {tTz/2 — b[hi + 62^2) and 



Hh 



fif^b ( hiTi- + b\Tr+ ) , if A > (28a) 



Hl^, = hil], ( bin- + b2TT+ ) , if A < . (28b) 



Here, 



f^b=5V(l-M)/(l + A^) 



|N+1> 



|N> 



|N-1> 






UJ2 



+ 



+ 



+ 



FIG. 2: Left ladder: energy of photon states |A'') (with 
A'^ S> 1) of the semiclassical field at frequency uj^. Mid- 
dle ladder: corresponding energies of the doublets of states 
{\g, N), je, N — 1)}, where A = ljl ~ 1^0 is the corresponding 
splitting in energy. Right ladder: energy of the semiclassical 
dressed states |±), Eq. (23), with energy splitting d, Eq. (22l. 
A transition |-|-) — )■ |— ) is accompanied by absorption (emis- 
sion) of a photon of frequency uji (1^2) from (into) the corre- 
sponding cavity mode. 



with = arctanh/i, while the value of fi is determined 



by the classical field parameters. 



/i = ta.n^e if |tan6l| < 1 
fi = {teine)-^ if Itane*! > 1. 



(29) 



We now can see that the interaction of the cavity modes 
with a beam of atoms, each initially prepared in the 
state |-|-) (|— )) and undergoing the dynamics governed 
by Hamiltonian H^ = Hq + H^^^^ {H2 ^ Hq + Ant)' '^iH 
give rise to an effective dynamics, whose steady state is 
the pure state lip) of the cavity modes. The Hamiltonian 
governing the dynamics will be Hi or H2 depending on 
the sign of the detuning A, which can be controlled by 
appropriately shifting the atomic transition frequency. 



B. Effective dynamics: reaching the EPR state 

We discuss now how the cavity modes can be prepared 
in the two-mode squeezed state asymptotically. This is 
achieved by an effective dissipation!^ process in the b- 
basis, implemented in a two-step procedure sketched in 
Fig. [3j The first step consists in letting the atomic beam 
interact with the resonator with each atom prepared in 
the state |-|-) and the detuning of the classical field set to 
the value A = Aq > 0. Lrsid e the cavity each atom un- 

such that at the end 



dergoes the dynamics of Eq. ( 28a 



of the interaction, on average, excitations have been re- 
moved from mode 61. In the second step, the atoms are 
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l+> 




J ^ — |g> 



C 0||| 



-|e> 



■|g> 



FIG. 3: Schematic representation of the interaction processes 
needed, in order to prepare the cavity modes in a EPR state. 
In step 1, the atoms are prepared in the state |+) before 
crossing the resonator. Inside the resonator they are driven 
by a classical field, and the atomic frequency is shifted with 
respect to lul according to the energy level scheme displayed 
at the bottom of the figure. In step 2, the initial state is |— ) 
and the atomic transition frequency is shifted, such that the 
detuning with the field has opposite sign. 



prepared in state |— ) and the detuning of the classical 
field is set to the value A — — Aq. Inside the resonator 



each atom undergoes the dynamics given by Eq. (28b), 



such that at the end of the interaction, on average, exci- 
tations have been removed from mode 62- 

In the weak-coupling regime, the equation of motion of 
the cavity field modes is given by the master equation for 
the density matrix pt, which during each step j {j — 1, 2) 
reads 



dpt 
dt 

where 



stop j 



('5 



hjpt - 2bjptb 



2_2 



7 = fat^^b T 



(30) 



(31) 



T is the interaction time, and rat is the atomic arrival 
rate. Correspondingly, during each step the average num- 
ber of Bogoliubov excitations is exponentially damped 



(6j6j)oexp(-7i), and 



according to the equation {Vjbj)t 
vanishes at times ^^1/7 (see also Sec. II C ) . In terms of 
the original field modes, this procedure implies that the 
atoms pump in phase only the two-mode squeezed state. 
Asymptotically, after the implementation of the two steps 
for a sufficiently long time (such that the cavity field is 
still stable over this time) the field state approaches the 
state 

Poo = |0, 0), (0, 0| - S\t^) |0, 0), (0, 0| S{r^), (32) 

which is a two-mode squeezed state, and whose degree 
of squeezing is solely determined by the ratio |A/ri|. 



This state is reached independently of the initial state of 
the cavity modes, provided that each step is implemented 
for a sufficiently long time T. 

An analogous dynamics can be implemented in the 
strong-coupling regime, i.e., when rfib ^ 1- Here, one 
aims at creating a trapping state in the Bogoliubov ba- 
sis for each Bogoliubov oscillator in a two-step procedure, 
where the initial state of the atomic beam and the atomic 
parameter are changed as in the procedure outlined for 
the weak-coupling case. In this case, the creation of trap- 
ping states in the Bogoliubov basis, with numbers rii 
and 712 for each Bogoliubov mode, corresponds to entan- 
gled states of the cavity modes of the form S{^)\ni,n2) ■ 
Choosing a broad velocity distribution, one finds that 
the steady state is a vacuum two-mode squeezed state, 
5(010,0). 



C. Experimental parameters 

The proposal we discussed so far is based on a stochas- 
tic dynamics, where the average action of each atom leads 
the cavity modes to a stationary, EPR-entangled state. 
The scheme does not require atomic detection, nor con- 
trol of the number of atoms, nor of the interaction times 
(atomic velocities). A possible experimental setup of this 
scheme is sketched in Fig. [T|a). Prior to the interaction 
region, the atoms are prepared in a coherent superposi- 
tion of two Rydberg states \g) and |e) connected by a 
dipole transition. Inside the resonator a classical field 
saturates the dipole transition, thereby pumping on res- 
onance the two nondegenerate modes of the resonator, as 
shown in Fig. [2] 

While the interaction time between cavity and each 
atom needs not to be controlled, on the other hand the 
dynamics between atom and cavity is here assumed to be 
Hamiltonian, and characterized by a two-level transition 
which can be tuned on resonance with the cavity modes 
by means of an external field. Let us now discuss these as- 
sumptions individually. Selecting a two-level transition, 
here denoted by the electronic states \g) — )■ |e) imposes 
constrains on the field polarizations, such that they all 
couple with the dipole transition, while coupling to other 
states is avoided. In an open-cavity geometry |22| . this 
can be achieved by means of an electric potential between 
the two mirrors, which removes through Stark shifting 
the degeneracy of circular Rydberg states, and using cir- 
cular polarizations for both the cavity modes and the 
pump field. The two-step procedure needs a change in 
the transition frequency of the two-level atom, which can 
be achieved by an external static field. 

The assumption of Hamiltonian dynamics between in- 
dividual atoms and resonator relies on the fact that the 
atom must not decay during the interaction with the cav- 
ity modes, and dissipation of the cavity field should be 
negligible during the experiment. Experiments with mi- 
crowave resonators |23] are characterized by inter- 
action times of the order of tens of ps, which warrant 
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negligible spontaneous decay, typically of the order of 
tens of ms. The requirement that the cavity does not 
decay over the duration of the experiment Ttot — 2T is 
instead more delicate, as it requires that the resonator 
is stable over the total interaction time with the atomic 
beam. Moreover, the time T must be sufficiently smaller 
than the coherence time of the driving fields, so that the 
amplitude f2 remains constant. 

Being the dynamics stochastic, the desired EPR-state 
is reached asymptotically. For this purpose, we estimate 
the time needed for reaching a given fidelity in the prepa- 
ration of the desired state. We first focus on the weak- 
coupling limit. Using Eq. ([9| we find that the fidelity of 
the protocol has the form 



1 



1 



(33) 



where we have assumed that each step is performed on 
the time interval T, such that the total duration of the 
protocol Ttot = 2T. Clearly, the initial state affects the 
time scale required for reaching the desired fidelity. In 
particular, when the cavity modes are initially in the vac- 
uum state one has {bjbj)o — —■ ^o- When they 

are in a thermal state with nth thermal photons in each 
mode, instead, thus (bjbj)o = tiq^ — uq + nth + '^n-thno- 
Figure [4] displays the estimated total experimental times 
and corresponding average number of photons per mode 
at steady state as a function of fi, where no = — 
when the cavity modes are in the vacuum state at i = 0. 
For the degree of squeezing w 2.1 ~ 0.97), lead- 
ing to an average number of 16 photons per mode at 
steady state, and = 0.01, corresponding to a fidelity 
0.98, then one has Ttot 36 ms in case of an ini- 
tially empty cavity (Ttot 43 ms for 0.7 thermal pho- 
tons). Resonators stable over this time are available in 
present experiments [28]. From these results we also see 
that fluctuations in the coupling with the driving field, 
dfl, are negligible with current microwave sources. 

We now consider the implementation of the protocol 
in the strong-coupling limit. In this case the fidelity for 
preparing the cavity modes in the two-mode squeezed 
state takes the form 



1 



(34) 



where we have used Eq. ( 15 1, assuming a broad distribu- 



tion of atomic velocity and an initially empty cavity. 

In order to compare the efficiency of the two proce- 
dures, we now evaluate the time required to reach a 
desired fidelity in the weak and in the strong-coupling 
limit. We first observe that the system accesses the 
strong-coupling regime by increasing the Rabi angle 4>. 
In the case here considered, (j) can be changed by chang- 
ing the interaction time while keeping g, and thus ilb 
fixed. Note that increasing the interaction time implies 
that the atomic arrival rate rat needs to be adjusted to 
avoid simultaneous presence of more than one atom in the 
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FIG. 4; (i) Solid line: Total experimetal time Ttot = 2T (in 
seconds) and (ii) Dashed line: average photon number per 
mode at the end of the protocol as a function of /i (which 
is controlled by the intensity of the classical field and the 
detuning according to Eq. (29) and (24)). The parameters 
we used are efi'ective coupling g — 125 kHz, interaction time 
f — 12.5/is, atomic arrival rate 11200 atoms/sec. 



cavity. Let us therefore take the average number of atoms 
in the interval of time t so that e = ratT ^ 1. In the 



weak-coupling limit the fidelity in Eq. ( 33 1 , shows that it 
is favorable to increase f even at the expense of a slower 
rate Tat (see Eq. (31 1). In the strong-coupling limit, from 
Eq. (34) we find that increasing f further gives slower 
convergence to the desired fidelity, hence slower proto- 
cols. Figure [5] displays the time required for reaching a 
desired fidelity — 0.99 as a function of f. One finds an 
optimal f, for which the protocol is fastest, in between 
the two regimes. 



IV. ENTANGLING TWO DISTANT CAVITIES 
WITH AN ATOMIC RESERVOIR 

Thus far we have considered two modes of the same mi- 
crowave cavity. However, as EPR entangled radiation is a 
possible resource for quantum communication with con- 
tinuous variables, e.g. for quantum teleportation [48], it 
would be desirable to create two entangled modes belong- 
ing to two different, spatially-separated cavities. In this 
section we extend the concepts discussed in Section III 
and present a scheme for entangling the modes of two 
spatially separated resonators using an atomic beam. 

We assume two open resonators, which are crossed by 
a beam of atoms propagating along the z axis and whose 
relevant modes are at frequency wi and uj2, with anni- 
hilation and creation operators ai,a\ and a2,a|, respec- 
tively. The resonators spatial mode functions along the 
z axis are Gaussians centered at zi and Z2, respectively, 
such that the distance |zi — Z2I — I? is much larger than 
the modes waist w, namely, D ^ w and the fields mode 
functions have no spatial overlap. The atomic transition 
is quasi-resonant with a mode at frequency uji in the first 
resonator, and at frequency uj2 with the second resonator. 
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FIG. 5: Logarithmic plot of time T required in each step to 
reach the EPR state as a function of average interaction time 
r (in units of f^,^^). The EPR state here chosen is charac- 
terized hy fi — .95 (no ~ 9.3), and we require that the cor- 
responding fidelity at Ttot = 2T is = 0.99. The blue lines 
represent analytical expressions for fidelities J^vc, Eq. and 
J-sc, Eq. (151. The red line has been obtained from numeri- 



cal solution of Eq. ( 10 1, where coefficients Bn were calculated 
taking a — 0.05<j!>o in the Gaussian distribution function of 
Rabi angles (j)o — fQt- The most efficient regime, in terms 
of fastest protocols for a given fidelity, lies between weak and 
strong coupling. The plot shows slower convergence at Rabi 
angles corresponding to <j)o — tt/v'j (j positive integer) due 
to trapping states. 



we will show, the atomic velocity can change without 
affecting the efficiency of the protocol. 

For the following treatment we assume that the mode 
waist is much smaller than the wavelength, as it is of- 
ten the case, so that the maser field has a well defined 
phase over the interaction region. In particular, we de- 
note by -ijji = kzi and V'2 = the phases at each res- 
onator. The atom is continuously driven by the maser 
field, and the phase of the dressed states is here assumed 
to follow adiabatically the phase of the field as the atom 
moves. For simplicity, in what follows we assume the 
phases V'l = "02 = 0. Within these approximations, we 
write the effective dynamics of the individual atom in- 
teracting sequentially with the cavity modes in terms of 
semiclassical dressed states of the maser field, given in 
Eqs. (23). Setting the detuning between maser and cav- 
ity modes lo^ ^ loi — d and — uj2 ~ —d, with d the 
frequency splitting between the dressed states given in 
Eq. (22), we obtain the corresponding effective Hamilto- 



nian, which reads HcS j = Hq 



Hj with 



Ho 
Hi 

H2 



M 



— a[ai + 0^02 



-figi sin^ 0fi{z) 



%2 cos 9/2(2) dl^TT +0,27: 



d\f[^ + fllTT 



(36a) 

(36b) 
(36c) 



In addition, the atoms are driven by a maser which prop- 
agates almost parallel to the z-axis and which has wave 
vector fc, frequency lul and intensity fl (we neglect any 
spatial gradient and assume that the maser intensity is 
uniform along z). 

The Hamiltonian describing the coherent interaction 
of one individual atom of the beam with both cavities 
reads 

H{t) = Hujoa^a (35) 
+ ^ l^fiwAaiflA + fi9\f\{z{t)){dxcr'^ + d\cr) , 

A 

where g\ is the strength of the coupling between cav- 
ity mode and dipolar transition, while the spatial mode 
function takes the form 

fx(z) = exp(-(z - za)V2u;')/^W, 

with w the mode waist. The atomic center of mass is 
located at the time-changing position z{t) — z(0) + vt, 
where v is the velocity of the atom. The assumption 
of classical center-of-mass motion is justified by the pa- 
rameters of the typical experimental situation, where the 
atoms exit an oven and the velocity selection brings to 
distributions still within the classical regime. In addi- 
tion, the assumption of uniform motion is not necessary, 
but convenient for the theoretical treatment. Indeed, as 



where the angle 9 is given in Eq. (24). The space- 
dependence of the cavity spatial mode functions can be 
substituted by the mean value of the function over the 
interaction region, which is constant over an interval of 
time, such that the total pulse area is preserved [16 . 



Hamiltonian Eq. (36) has the same form as Eq. (27) in 
Sec. Ill with the only difference that the atom interacts 
with each mode in well-separated time intervals, due to 
the distant location of the cavities. We will now show 
that the combined state of the two distant cavity modes 
can be pulled into a two-mode squeezed state as a result 
of the interaction with an atomic beam, of which the 
number of atoms and the individual interaction times 
are known only statistically. 

The line of reasoning extends the protocol presented 
in Sec. |III[ In this case, however, one must consider that 
the interaction of the atom with each mode is sequen- 
tial. For this purpose we introduce the evolution opera- 
tor Uj^r) for the dynamics of the atom interacting with 
the resonator j over the interval of time t (here taken 
in interaction picture with respect to Hamiltonian Hq in 
Eq. (36a)) and we evaluate the total density matrix pj^ of 
atom and resonator after the atom has interacted with 
both cavities in interaction picture, assuming that the 
atom interacts sequentially first with mode 1 and then 
with mode 2. Denoting by At the interval of time in 
which this occurs, the density matrix after the interac- 
tion reads 



nt + At) = U2{T)Ul{T)p{t)Ul{T)Ul{T) 



(37) 
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In order to obtain the two-mode squeezing correlations, 
we need also processes in which the temporal sequence 
of photon absorption and emission is reversed, namely, 
processes of the sort 



(38) 



and which require the presence of a second current of 
atoms propagating in the opposite direction, such that 
each atom first interacts with resonator 2 and then with 
resonator 1, see Fig. |6| After imposing this condition, 
one derives an effective dissipative dynamics fo r th e Bo- 
goliubov modes hi and 62, defined as in Eq. (18), and 



which is valid over a time step Ai such that there is at 
most one atom inside the resonator at a time. Provided 
that the atoms are initially prepared in the state |+) 
(|— )) and the detuning is such that A = Aq (A = — Aq), 
the master equation for the field density matrix in the 
coarse-grained time-scale is 



dt 
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b]b,p 



2bjpb 



(39) 



where Ap^ = p^{t + At) — p{t) (same for ^), and rat 
is the atom pump rate, which is assumed to be the same 
in both directions. Hence, mode bi (62) is exponentially 
damped according to an equation of the form given in 
Eq. (30). Preparation of the two resonators in a two- 



mode squeezed state then is achieved, provided that both 
dynamics take place, by changing detuning and state 
preparation after the first step has been implemented af- 
ter a sufficiently long time. 

As in the protocol for entangling two modes of the same 
resonator, this scheme does not require detection of final 
atomic states, nor control of the atomic velocities. More- 
over, the velocity of the atoms can change during prop- 
agation, without affecting the efficiency of the protocol. 
This result has been derived considering Hamiltonian dy- 
namics between cavity and individual atoms. 

Let us now discuss the limitations to this proposal. 
First, only one atom at a time must be present inside the 
cavities. Moreover, the atoms must not decay before they 
have interacted with both resonators. Atomic lifetime of 
the order of tens of milliseconds and atomic velocities 
of about 400 m/s require that the distance between the 
cavities is no larger than a few meters. Moreover, the 
resonators must be stable over the whole run of the ex- 
periment. 

Experimental implementation of this proposal is quite 
challenging, since one must have a geometry with two 
counterpropagating atomic beams, thereby avoiding col- 
lisions between the atoms. As an alternative, one could 
think of a ring, as realized for instance in ion-storage se- 
tups |51) in which the atomic beam is confined and which 
crosses two resonators placed at two different points of 
the ring. Another possibility is to implement an atomic 
fountain, where atoms traveling upwards and downwards 



l+> 




-l+> 




+> 
i-> 



FIG. 6; Top: sketch of the setup for entangling the modes of 
two distant resonators. Two atomic beams propagate in both 
directions and cross the resonators, interacting sequentially 
with each of them. Here one of the two required steps of 
the procedure is shown, where the atoms are prepared in the 
state |-|-). Bottom: energy levels which are coupled inside 
each resonator. The "ball" represents the atomic occupation. 
An atom traveling from right to left first emits a photon UJ2 
into the second resonator, then emits a photon wi into the 
first one, only if it has previously emitted into resonator 2, 
otherwise it may absorb a photon uii . An atoms traveling from 
left to right absorbs/emits photons in the reversed sequence. 



interact with two vertically arranged cavities above the 
fountain [52], or a optical conveyor belt [53], where the 
atoms are transported back and forth between the res- 
onators. These mechanisms have been employed so far 
in the optical regime, where the lifetime of the resonator 
modes is limited to tens of microseconds, in which case 
the implementation of this protocol is not straightfor- 
ward. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have extensively characterized the 
properties of an atomic beam as a reservoir for the 
modes of the electromagnetic field inside a resonator, for 
the purpose of creating entangled states of the cavity 
modes. The atoms can mediate the interaction between 
the modes of the same resonator or of distant cavities, es- 
tablishing a dynamics whose steady state is an Einstein- 
Poldosky-Rosen entangled state. As opposed to previous 
proposals, see for instance [36l |37], the atoms do not 
need to be initially correlated nor their number has to 
be controlled. Control on atomic velocity (interaction 
time) and atomic detection are not required. The de- 
gree of entanglement is controlled by an external maser 
field, which drives the atoms and tailors their interaction 
with the cavity modes. In this respect, the proposals 
discussed in this article are instances of quantum reser- 
voir engineering. Statistical properties of the cavity field 
can be evaluated by measuring the internal states of the 
emerging atoms [13]. Its state can also be determined 
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by reconstructing the corresponding Wigner function, by 
suitably generalizing the schemes proposed in |50j . 

The experimental setup, where these protocols could 
be implemented, is typical of microwave cavity quantum 
electrodynamics. Our proposal extends to cavity QED 
the technique of quantum reservoir engineering, origi- 
nally applied to trapped ions. It offers a convenient way 
of implementing quantum mechanical dynamics and state 
preparation. Besides constituting a robust procedure for 
generating nonclassical states of the electromagnetic field 
in cavities in a steady-state regime, it might be useful for 
quantum networking with continuous variables in the mi- 
crowave regime. 
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